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(54) Apparatus and method for removing scatter from an x-ray image 



(57) An apparatus and method for eliminating scat- 
ter effects in x-ray imaging using two-dimensional de- 
tector arrays is disclosed. The apparatus consits of an 
x-ray source (14), a front two-dimensional x-ray detector 
(16), a beam selector (18). and a rear two-dimensional 
x-ray detector (26). A substantially improved dual-ener- 
gy method is presented. Based on this method, the pri- 



mary image of the front detector (16) is calculated from 
a pair of primary dual-energy images of the rear detector 
(26). Subsequently, the low resolution scatter image of 
the front detector (16) is determined and the high reso- 
lution scatter image of the front detector (16) is calcu- 
lated. Embodiments of this invention provide a method 
for elimination of scatter and a method for dual energy 
x-ray imaging. 
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Description 



This invention relates generally to digital x-ray imaging and, more particularly, relates to methods and apparatuses 
for reducing scatter in two-dimensional x-ray imaging and two-dimensional dual-energy x-ray imaging. 

Recent advances in semiconductor industry have resulted in the ability of fabricating large-format two-dimensional 
integrated detector arrays for x-ray detection. These arrays have on the order of one million detector cells and provide 
instant acquisition of two-dimensional x-ray images with exceedingly high quality. 

Scatter, which results from those x-rays that strike objects and deflect in random directions, has been a difficult 
and on-going problem in x-ray imaging using two-dimensional detectors. For example, in projection chest radiography, 
scatter typically accounts for between approximately 30% and 50% of the total amount of x-rays detected. 

Two-dimensional detectors are exposed to wide-angle random scatter. Scatter tends to degrade image quality, 
and to block the way for quantitative imaging. For example, with scatter present, three-dimensional image reconstruc- 
tion cannot be implemented. 

Embodiments of the present invention use the dual-energy x-ray imaging method as a foundation for the method 
of eliminating scatter. However, the dual-energy x-ray imaging methods of prior art are not suitable for this purpose. 

AH the dual-energy x-ray data decomposition methods of the prior art have a common approach: linearization. 

It is a consensus that none of current linearization methods is suitable for use with million-pixel two-dimensional 
detector arrays. One object of the invention is to provide a substantially improved dual-energy x-ray imaging method 
suitable for two-dimensional detector arrays. 

Currently, there are two basic methods for reducing scatter in two-dimensional x-ray imaging. The first method 
uses an anti-scatter grid or improved grid devices to slightly relieve scatter effects on images. By increasing the air 
gap between the subject and the detector, the scatter can also be reduced, but at the same time, the image is blurred 
due to the geometric distance the x-rays have to travel. 

The second method to reduce scatter is to calculate theoretical estimates. Theoretical methods, including Monte 
Carlo simulation methods and analytical deconvolution methods, can only give crude predictions. 

According to a first aspect of the present invention, there is provided a two-dimensional x-ray imaging system for 
taking images of a subject, said system comprising: 

(a) an x-ray source adapted to emit x-rays for passage through said subject, said x-rays including primary x-rays 
having their direction of travel unaltered by interaction with said subject and said x-rays including scatter x-rays 
having their direction of travel altered by interaction with said subject; 

(b) x-ray detection means for detecting said x-rays in two dimensions; 

(c) said subject being located between said x-ray source and said x-ray detection means; 

(d) said detection means including high-resolution detection means for detecting said primary x-rays and said 
scatter x-rays simultaneously; 

(e) beam selection means for blocking passage of substantially all of said scatter x-rays and permitting passage 
of a portion of said primary x-rays; 

(f) said detection means including low-resolution detection means for detecting only said x-rays passed by said 
beam selection means. 



According to a second aspect of the present invention, there is provided a method for taking a twondimensional x- 
ray image of a subject using a two-dimensional x-ray imaging system, said subject being composed substantially of 
two materials that interact differently with x-rays, said imaging system including an x-ray source, beam selection means, 
and a two-dimensional x-ray detector assembly having a plurality of detection locations, said subject being between 
said x-ray source and said detector assembly, said x-ray source being adapted to emit x-rays for passage through said 
subject, said x-rays including primary x-rays having their directbn of travel unaltered by interaction with said subject, 
said x-rays including scatter x-rays having their direction of travel altered by interaction with said subject, said beam 
selection means having a blocking position wherein passage of substantially all of said scatter x-rays to said detector 
assembly is blocked and passage of said primary x-rays to said detector assembly is permitted, said beam selection 
means having a transmissive position wherein passage of substantially all of said primary x-rays and said scatter x- 
rays to said detector assembly is permitted, and selected detection locations being those of said detection locations 
receiving only said primary x-rays when said beam selection means is in said blocking position, said method comprising: 

(a) illuminating said subject with said x-rays; 

(b) moving said beam selection means to said transmissive position; 

(c) acquiring a raw high-resolution image l m from said detection locations and processing said image l^ to normalize 
it and to subtract dark signals, yielding a high-resolution image D m composed of said primary x-rays and said 
scatter x-rays; 
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(d) producing, from said image D^, a low- resolution image D n representing said selected detection locations; 

(e) moving said beam selection means to said blocking position; 

(f) acquiring a raw low- resolution image l rt from said selected detection locations and processing said image l ri to 
normalize it and to subtract dark signals, yielding a low-resolution image D rt composed ol substantially only said 

5 primary x-rays; 

(g) calculating a low-resolution primary x-ray image D fP , by multiplying said image D ri by a predetermined constant 
that is independent of said subject; 

(h) calculating a low-resolution scatter x-ray image Dfsi by subtracting said image Dfp, from said image D„; 

(i) calculating a high-resolution scatter x-ray image Dfs h by extending said low-resolution scatter x-ray image 
10 to all of said detection locations through interpolation; and 

(j) calculating a high-resolution primary x-ray image Dfp h by subtracting said high-resolution scatter x-ray image 
D fSh from said high-resolution image D^; 

(k) whereby said image Dfp h is a high-resolution, two-dimensional primary x-ray image of said subject after said 
scatter x-rays have been substantially eliminated, said image Dfp h having a resolution substantially equal to the 
is highest resolution available from said detector assembly. 

According to a third aspect of the present invention, there is provided a method for taking a two-dimensional x-ray 
image of a subject using a two-dimensional x-ray imaging system, said subject being composed substantially of two 
materials that interact differently with x-rays, said imaging system including, in physical sequence from front to back, 

20 an x-ray source, a front two-dimensional x-ray detector assembly having a plurality of front detection locations, beam 
selection means, and a rear two-dimensional x-ray detector assembly having a plurality of selected rear detection 
locations, said subject being between said x-ray source and said front detector assembly, said x-ray source being 
adapted to emit x-rays for passage through said subject, said x-rays having an energy spectrum that can be approxi- 
mated as having a single energy E 0 , said x-rays including primary x-rays having their direction of travel unaltered by 

25 interaction with said subject and said x-rays including scatter x-rays having their direction of travel altered by interaction 
with said subject, said beam selection means blocking the passage of substantially all of said scatter x-rays to said 
selected rear detection locations and permitting the passage of said primary x-rays to said selected rear detection 
locations, said selected rear detection locations receiving only primary x-rays, and selected front detection locations 
being those of said front detection locations intersected by x-ray projection lines extending from said x-ray source to 

30 ..said selected rear detection locations, said method comprising: 

(a) .illuminating said subject with said x-rays; 

(b) acquiring a raw high-resolution image l ft from said front detection locations and processing said image l m to 
normalize it and to subtract dark signals, yielding a high-resolution image D m composed of primary x-rays and 

35 scatter x-rays; 

(c) producing, from said image D mt a low- resolution image D fl representing said selected front detection locations; 

(d) acquiring a raw low-resolution image l rt from said selected rear detection locations and processing said image 
l rt to normalize it and to subtract dark signals, yielding a low-resolution image D rt composed of substantially only 
said primary x-rays; 

40 (e) calculating a low- resolution primary x-ray image D fP , at said front detector assembly by multiplying said image 

D rt by a predetermined constant that is independent of said subject; 

(f) calculating a low-resolution scatter x-ray image D fe , at said front detector assembly by subtracting said image 
D fP | from said image D f) ; 

(g) calculating a high-resolution scatter x-ray image D^ h by extending said low-resolution scatter x-ray image D^, 
45 to all of said front detection locations through interpolation; and 

(h) calculating a high-resolution primary x-ray image D fPh by subtracting said high-resolution scatter x-ray image 
D fSh from said high-resolution image D^; 

(i) whereby said image Dfp h is a high-resolution, two-dimensional primary x-ray image of said subject after said 
scatter x-rays have been substantially eliminated, said image Dfp h having a resolution substantially equal to the 

so highest resolution available from said front detector assembly. 

According to a fourth aspect of the present invention, there is provided a method for taking a two-dimensional x- 
ray image of a subject using a two-dimensional x-ray imaging system, said subject being composed substantially of 
two materials that interact differently with x-rays, said imaging system including, in physical sequence from front to 
55 back : an x-ray source, a front two-dimensional x-ray detector assembly having a plurality of front detection locations, 
beam selection means, and a rear two-dimensional x-ray detector assembly having a plurality of selected rear detection 
locations, said subject being between said x-ray source and said front detector assembly, said x-ray source being 
adapted to emit x-rays at two different average energy levels for passage through said subject, said x-rays including 
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primary x-rays having their direction of travel unaltered by interaction with said subject and said x-rays including scatter 
x-rays having their direction of travel altered by interaction with said subject, said beam selection means blocking the 
passage of substantially all of said scatter x-rays to said selected rear detection locations and permitting the passage 
of said primary x-rays to said selected rear detection locations, said selected rear detection locations receiving only 
s primary x-rays, and selected front detection locations being those of said front detectbn locations intersected by x-ray 
projection lines extending from said x-ray source to said selected rear detection locations, said method comprising: 

(a) illuminating said subject with x-rays of a first average energy level H; 

(b) acquiring a raw high-resolution image l^ from said front detection locations and processing said image I^k 
io to normalize it and to subtract dark signals, yielding an image D fHh which is composed of primary x-rays and scatter 

x-rays; 

(c) producing, from said image D mh , a low-resolution image D^,, representing said selected front detection loca- 
tions; 

(d) acquiring a raw low- resolution image 1^, from said selected rear detection locations and processing said image 
15 'rHi to normalize it and to subtract dark signals, yielding an image D^, composed of substantially only said primary 

x-rays; 

(e) illuminating said subject with x-rays of a second average energy level L; 

(f) acquiring a raw high-resolution image 1^ from said front detection locations and processing said image 1^ to 
normalize it and to subtract dark signals, yielding an image which is composed of primary x-rays and scatter 

zo x-rays; 

(g) producing, from said image D^, a low-resolution image Dflj representing said selected front detection locations; 

(h) acquiring a raw low-resolution image 1^, from said selected rear detection locations and processing said image 
IrL, to normalize it and to subtract dark signals, yielding an image D^, composed of substantially only said primary 
x-rays; 

2S (j) calculating a pair of low-resolution primary x-ray images D^p, and D^p, at said front detector assembly from a 

low-resolution dual-energy primary x-ray imaging pair composed of said image D^, and said image D,^,; 
(j) calculating a low- resolution scatter x-ray image D^ SI at said front detector assembly by subtracting said image 
D^p, from said image and calculating a low- resolution scatter x-ray image D^q, at said front detector assembly 
by subtracting said image D^p, from said image D rHI ; 

30 (k) calculating a high-resolution scatter x-ray image D^ Sh by extending said low-resolution scatter x-ray image 

DfLQ, to all of said front detection locations through interpolation and calculating a high-resolution scatter x-ray 
image DfH S h by extending said low-resolution scatter x-ray image D fHS i to all of said front detectbn locations through 
interpolation; and 

(I) calculating a high -resolution primary x-ray image at said front detector assembly by subtracting said image 
35 D fHSh trom saj d image D fHh and calculating a high-resolution primary x-ray image at said front detector 

assembly by subtracting said image D fLSh from said image D^; 

(m) whereby said images D fHPh and D,^ form a high-resolution, two-dimensional, dual-energy primary x-ray 
image pair D mPh , of said subject after said scatter x-rays have been substantially eliminated, said image 
pair having a resolution substantially equal to the highest resolution available from said front detector assembly. 

40 

According to a fifth aspect of the present invention, there is provided a method for taking a two-dimensional x-ray 
image of a subject using a two-dimensional x-ray imaging system, said subject being composed substantially of two 
materials that interact differently with x-rays, said imaging system including, in physical sequence from front to back, 
an x-ray source, a front two-dimensional x-ray detector assembly having a plurality of detectbn locations, beam se- 
4S lection means, and a rear detector assembly assembly, said subject being between said x-ray source and said front 
detector assembly, said x-ray source being adapted to emit x-rays for passage through said subject, said x-rays in- 
cluding primary x-rays having their direction of travel unaltered by interaction with said subject and said x-rays including 
scatter x-rays having their direction of travel altered by interaction with said subject, said rear detector assembly in- 
cluding a rear low-energy two-dimensional x-ray detector having a plurality of selected rear low-energy detection lo- 
cations, a rear x-ray energy spectral filter, and a rear high-energy two-dimensional x-ray detector having a number of 
selected rear high-energy detection locations substantially equal to said plurality of selected rear low-energy detection 
locations, said beam selection means blocking the passage of substantially all of said scatter x-rays to said selected 
rear bw-energy detection locations and said selected rear high-energy detection locations and permitting the passage 
of said primary x-rays to said selected rear low-energy detection locations and said selected rear high-energy detection 
locations, said selected rear low-energy detection locations and said selected rear high-energy detection locations 
receiving only primary x-rays, and selected front detection locations being those of said front detection locations inter- 
sected by x-ray projection lines extending from said x-ray source to said selected rear low-energy detection locations 
and said selected rear high-energy detection locations, said method comprising: 
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(a) illuminating said subject with said x-rays; 

(b) acquiring a raw high -resolution image l^ from said front detection locations and processing said image l m to 
normalize it and to subtract dark signals, yielding a high-resolution image D m , which is composed of said primary 
x-rays and said scatter x-rays; 

5 (c) producing, from said image D m , a low-resolution image D fl representing said selected front detection locations; 

(d) acquiring a raw low- resolution image l^, from said selected rear high-energy detection locations and processing 
said image l^, to normalize it and to subtract dark signals, yielding a low-resolution image D^i composed of 
substantially only said primary x-rays; 

(e) acquiring a raw low-resolution image l^j from said selected rear low-energy detection locations and processing 
10 said image l^, to normalize it and to subtract dark signals, yielding a low-resolution image composed of sub- 
stantially only said primary x-rays; 

(f ) calculating a low-resolution primary x-ray image D m at said front detector assembly from a low-resolution dual- 
energy primary x-ray imaging pair composed of said image D^, and said image D^; 

(g) calculating a low-resolution scatter x-ray image Dfs, at said front detector assembly by subtracting said low- 
1S resolution primary x-ray image Dfp, from said low-resolution image D fl ; 

(h) calculating a high-resolution scatter x-ray image by extending said low- resolution scatter x-ray image D^, 
to ail of said front detection locations through interpolation; and 

(i) calculating a high-resolution primary x-ray image Dfp h by subtracting said high-resolution scatter x-ray image 
Dfs h from said high-resolution image D m ; 

20 (j) whereby said high-resolution primary x-ray image D fPh is a two-dimensional image of said subject after said 

scatter x-rays have been substantially eliminated, said image Dfp h having a resolution substantially equal to the 
highest resolution available from said front detector assembly. 

According to a sixth aspect of the present invention, there is provided a method for taking a two-dimensional x- 

25 . ray image of a subject using a two-dimensional x-ray imaging system, said subject being composed substantially of 
two materials that interact differently with x-rays, said imaging system including, in physical sequence from front to 
„ back, an x-ray source, a front detector assembly, beam selection means, and a rear detector assembly, said subject 
being between said x-ray source and said front detector assembly, said x-ray source being adapted to emit x-rays for 
passage through said subject, said x-rays including primary x-rays having their direction of travel unaltered by inter- 

30 , action with said subject and said x-rays including scatter x-rays having their direction of travel altered by interaction 
with said subject, said front detector assembly including a front low-energy two-dimensional x-ray detector having a 
.plurality of front low-energy detection locations, a front x-ray energy spectral filter, and a front high-energy two-dimen- 
sional x-ray detector having a number of front high-energy detection locations equal to said plurality of front low-energy 
detection locations, said rear detector assembly including a rear low-energy two-dimensional x-ray detector having a 

35 plurality of selected rear low-energy detection locations, a rear x-ray energy spectral filter, and a rear high-energy two- 
dimensional x-ray detector having a number of selected rear high-energy detection locations equal to said plurality of 
selected rear low-energy detection locations, said beam selection means blocking the passage of substantially all of 
said scatter x-rays to said selected rear low-energy detection locations and said selected rear high-energy detection 
locations and permitting the passage of said primary x-rays to said selected rear low-energy detection locations and 

40 said selected rear high-energy detection locations, said selected rear low-energy detection locations and said selected 
rear high-energy detection locations receiving only primary x-rays, selected front high-energy detection locations being 
those of said front high-energy detection locations intersected by x-ray projection lines extending from said x-ray source 
to said selected rear low-energy detection locations and said selected rear high-energy detection locations, selected 
front low-energy detection locations being those of said front low-energy detection locations intersected by said x-ray 

45 projection lines, said method comprising: 

(a) illuminating said subject with said x-rays; 

(b) acquiring a raw high-resolution image l^ from said front low-energy detection locations and processing said 
image l^ to normalize it and to subtract dark signals, yielding a high-resolution image Dfl_ h , which is composed 

50 of said primary x-rays and said scatter x-rays; 

(c) producing, from said image Dfl_ h , a low-resolution image Dfy representing said selected front low-energy de- 
tection locations; 

(d) acquiring a raw high-resolution image l^ from said front high-energy detection locations and processing said 
image \ mu to normalize it and to subtract dark signals, yielding a high-resolution image D^, which is composed 

5 $ of said primary x-rays and said scatter x-rays; 

(e) producing, from said image D mh , a low- resolution image D^, representing said selected front high-energy 
detection locations; 

(f) acquiring a raw low-resolution image l^from said selected rear high-energy detection locations and processing 
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said image l rH , to normalize it and to subtract dark signals, yielding a low-resolution image D rHI composed of 
substantially only said primary x-rays; 

(g) acquiring a raw low-resolution image l^j from said selected rear low-energy detection locations and processing 
said image l^, to normalize it and to subtract dark signals, yielding a low-resolution image D^j composed of sub- 
stantially only said primary x-rays; 

(h) calculating a pair of low-resolution primary x-ray images D^p,; at said front high-energy detector and at 
said front low-energy detector from a low-resolution dual-energy primary x-ray imaging pair composed of said 
image D rH , and said image D,^,; 

(i) calculating a low-resolution scatter x-ray image D msi at said front high-energy detector by subtracting said 
image D^p, from said image D^j and calculating a low-resolution scatter x-ray image D^q, at said front low-energy 
detector by subtracting said image D fLP , from said image D^; 

(j) calculating a high-resolution scatter x-ray image by extending said low-resolution scatter x-ray image 
D fHSi t0 a!l of sai d front high-energy detection locations through interpolation and calculating a high-resolution 
scatter x-ray image D fLSh by extending said low-resolution scatter x-ray image Df|_ s , to all of said front low-energy 
detection locations through interpolation; and 

(k) calculating a high -resolution primary x-ray image D^ph by subtracting said image D fHSh from said image 
and calculating a high-resolution primary x-ray image D^p,, by subtracting said image D,^ from said image D^; 
(I) whereby said images D mPh and form a high-resolution, two-dimensional, dual-energy primary x-ray image 
pair D^ph, D^ph of said subject after said scatter x-rays have been substantially eliminated, said image pair having 
a resolution substantially equal to the highest resolution available from said front detector assembly. 

The dual-energy x-ray data decomposition method in this invention consists in directly solving the dual energy 
nonlinear equation system, it includes: (1) Constructing an explicit quantitative equation system Dh ~ D H (b,s) and D L 
= D L (b,s) for each detector according to the nonlinear dual-energy x-ray imaging fundamental equation system in its 
original form. (2) reconstructing an equation system b = b(D H ,D L ) and s = s(D H ,D L ) by numerically inverting the equation 
system of step 1.(3) Determining the desired values for b and s from the available data pair (Dh,D l ) by using the 
numerical equations of step 2, or determining the desired values for D H ,D L , from the available data pair (b,s) by using 
the numerical equations of step 1 . (4) Maintaining the accuracy at each step to be as high as real number analytical 
solutions can provide. 

Embodiments of the present invention advantageously comprise an x-ray detection system that has the following 
features: (1) using a beam selection means to physically separates the primary x-rays from the scatter x-rays at a 
number of selected locations; (2) using a front detector assembly detects the primary and scatter x-rays, using a rear 
detector assembly detects only the primary x-rays at selected locations without scatter. 

The method for eliminating scatter x-rays advantageously includes: (1) using a dual-energy method to determine 
the primary image of the front detector at the corresponding selected locations from the image data of the rear detector 
assembly; (2) using an interpolation method to extend the scatter component of the front detector image determined 
at a the selected locations to the entire image area of the front detector. Following this, a complete separation of scatter 
and primary images for the front detector at high spatial resolution is achieved. 

Thus, embodiments of the present invention seek to provide an apparatus and method for substantially eliminating 
the effects of scatter on two-dimensional x-ray detectors. 

Embodiments of the present invention may also seek to provide an apparatus and a method for scatter-free dual- 
energy x-ray imaging using two-dimensional detectors and to provide two high accuracy material composition images 
of a subject at the spatial resolution as high as a two-dimensional detector array can provide. 

For a better understanding of the present invention, and to show how it may be carried into effect, reference shall 
now be made by way of example to the accompany drawings, in which: 

Fig. 1 is a diagram of the basic hardware of an embodiment of the present invention; 

Fig. 2 is a curve describing a typical x-ray source energy spectrum used in embodiments of the present invention; 
Fig. 3 is a flow diagram of the basic dual-energy data decomposition method using the hardware of Fig. 1 ; 
Fig. 4 is a flow diagram of the basic method to eliminate scatter using the hardware of Fig. 1 ; 
Fig. 5 is a diagram of a first embodiment of the present invention; 

Fig. 6 is a flow diagram of the method of the first embodiment using the hardware of Fig. 5; 
Figs. 7a and 7b are diagrams of a special case of the first embodiment of the present invention; 
Fig. 8 is a flow diagram of the method of the special case using the hardware of Figs. 7a and 7b; 
Fig. 9 is a diagram of a second embodiment of the present invention; 

Fig. 10 is a flow diagram of the method of the second embodiment using the hardware of Fig. 9; 
Fig. 11 is a diagram of a third embodiment of the present invention; 

Fig. 12 is a flow diagram of the method of the third embodiment using the hardware of Fig. 11 ; and 
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Figs. 1 3a to 1 3d is a graphically representation of a method for inverting the nonlinear dual-energy equation sys- 



The basic apparatus includes five components: (1) an x-ray source, (2) a front two-dimensional x-ray detector 
assembly, (3) a beam selector, (4) a rear two-dimensional x-ray detector assembly, and (5) a computer. 

As shown in Fig. 1 , the subject under examination 1 2 is located between the x-ray source 1 4 and the front detector 
assembly 16. The x-ray source 14 emits x-rays 30, preferably with an energy in the range of from 10 keV to 500 keV, 
10 as shown in Fig. 2. A portion of the x-rays 32 passes through the subject 12 directly to the front detector assembly 16 
without a change in their direction of propagation. These x-rays 32 are called the primary x-rays and convey true 
information about the subject 1 2. The remainder of the x-rays 34 are randomly scattered as a result of interaction with 
the material of the subject 12. These x-rays 34 are called scatter and cause a distortion of the true information. 

The front detector in the front detector assembly 16 has a number of detector types. The first uses thin film amor- 
15 phous silicon as photodetection medium. An scintillation screen, is placed in close contact with the entire photosensitive 
area of the photodetector array. Preferably, the detector array has dimensions of 20 centimeters (cm) by 20 cm or 40 
cm by 40 cm for a single detector module. The cell size for this detector array is in the range ol from approximately 50 
|im by 50 um to approximately 1 mm by 1 mm. 

A second type of detector array uses an amorphous selenium film or selenium alloy film as the x-ray sensitive 
20 medium. Other typical two-dimensional detector arrays include but not limited to, charge-couple device (CCD) detec- 
tors, thin-film thallium -bromide-based detector arrays, avalanche silicon detector arrays, and phosphorstimulatable 
computed radiography screens. 

The front detector assembly 16 conveys an image that is a combination of the primary x-rays 32 and the scatter 
34, and is denoted by 



where D f denotes an image in the front detector assembly 1 6 and (x.y) denotes the two-dimensional Cartesian 
coordinates of a cell of the front detector assembly 16. D fPh (x,y) denotes the contribution from the primary x-rays 32 
and D fSh (x,y) denotes the contribution from the scatter 34. 

The beam selector is sandwiched between the front detector assembly 16 and the rear detector assembly 26. A 
preferred embodiment of the x-ray beam selector 18 is a quantity of x-ray-absorbent material having a large number 
of holes 20. The holes 20 are fabricated such that their axes are aligned with the travel direction of the primary x-rays 
32. As a result, the holes 20 permit all x-rays traveling along the axes of the holes 20 to pass through, while all x-rays 
traveling in directions deviating slightly from the hole axes are completely absorbed by the bulk material, of the beam 
selector 18. Because the holes 20 always have a finite size, a small portion of randomly scattered x-rays from the 
image subject 12 still can reach the rear detector assembly 26. However, as long as the hole size 20 is small and the 
thickness of the beam selector 18 is sufficiently large, this portion of the scatter 34 can be controlled to be negligibly 
small. Preferably, the holes 20 have a diameter that is in the range of approximately from 0.5 mm to 10 mm. Preferably 
there are as many holes as possible in the beam selector 1 8. However, the beam selector material must occupy enough 
area so that the amount of scatter 34 is reduced to an insignificant quantity. A compromise based on these factors 
results in a pitch that is preferably between 2 mm and 50 mm. Preferably, the x-ray source 14 is located between 20 
cm and 150 cm from the rear surface 24 of the beam selector 18. The invention holds equally true when the x-ray 
source is a point source or has a finite size. 

The term "selected location" is defined as the location on the rear detector assembly 28, or on a rear detector, 
where, due to the function of the beam selector, only primary x-rays are received, and the scatter x-rays are substantially 
blocked. The "selected projection line* is defined as a straight line connecting the x-ray source to a point in the "selected 
location". The rear detector cells at the selected locations have a fixed geometric relation with some of the front detector 
cells. This relation is established by drawing a selected projection line from the x-ray source 14 through the beam 
selector 18 to the selected location. This selected projection line intersects the rear detector surface at a rear detector 
cell at a coordinate (i,j), and intersects the front detector surface at a front detector cell at a coordinate (x(i),y(j)). Here 
W')>y(j)) denote the Cartesian coordinate (x,y) of the front detector cell in the front detector assembly 1 6 closest to the 
selected projection line. The relationship between (i,j) and (x(i),y(j)) is experimentally established for all of the holes 
20 of the beam selector 18 and stored in computer The images represented by the composite of the signals from the 
detector cells of only one detector on the selected projection lines are low-resolution images and are represented by 
the subscript lower-case 1. The images represented by the composite of the signals from all the front detector cells 
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D fh (x,y) = D fph (x > y) + D fSh (x,y) 



(1) 
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are high-resolution images and are represented by the subscript lower-case h. 

In connection with the material composition of the image subject, four quantities are defined. b(i,j) and s(i,j) are 
defined as the selected projection mass densities along the selected projection line b(x,y) and s(x,y) are defined 
as the projection mass densities along the projection line (x,y). The "projection mass density" is defined as the integrated 
s total mass of the image subject along the projection line per unit area. Because the projection mass density is not 
dependent on the size of detector cells, b(x(i),y(j)) = b(i,j) and s(x(i),y(j)) = s(i,j). 

The method for eliminating scatter, shown in the flow diagram of Fig. 4, consists of the following steps: (1 ) illumi- 
nating the subject with x-rays from the x-ray source 1 4; (2) acquiring a low-resolution image pair I^KU ) and 1^,0, j) from 
the rear detector assembly and processing it to normalize and to subtract dark signals, yielding a low-resolution image 
io pair D rH ,(i, j) and D rtJ (i, j) that are functions of the subject materials, where "acquiring an image" is defined as transferring, 
via electronic control circuits, the electrical signals induced by the x-ray illumination on each detector cell from a detector 
array to computer memory; (3) solving the image pair D rHI (i : j) and D^^ij) to determine the selected projection mass 
density of the subject materials b(i,j) and s(i,j) according to dual energy data decomposition methods; (4) acquiring a 
high-resolution image ^(x.y) from the front detector and processing it to normalize and to subtract dark signals, yielding 
a high-resolution image D^x.y), which is the sum of primary x-rays and scatter x-rays; (5) calculating the low-resolution 
primary image of the front detector D f p,(x(i),y(j)) at the detector cells on the selected projection lines from the mass 
projection densities b(i,j) and s(i,j); (6) subtracting DfprfxOJ.y 0)) from D^(x{\) y y{\)) to calculate the low-resolution scatter 
component D fS ,(x(i),y(j)) of the image D^xO^yfj)) at the detector cells on selected projection lines; (7) interpolating 
DfsK x 0)(x(j)) tor those front detector cells not on the selected projection lines : yielding the high-resolution scatter image 
D fSh( x >y); (8) subtracting the image D fSh (x,y) from D^x.y) to yield an image Dfp^x.y), which is a full two-dimensional 
image of the subject at the front detector after scatter x-rays have been substantially eliminated. 

Because the data decomposition method of the present invention is rigorous and accurate and the mathematical 
calculations can be performed with good precision, as long as the experimental image data acquired from the detectors 
16, 26 is highly accurate, the final result will be a highly accurate primary image of the subject 12. 
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In the first embodiment of the apparatus, shown in Fig. 5, the rear detector assembly 26 is constructed as a dual- 
energy x-ray imaging detector assembly. It has a low-energy two-dimensional detector 40, an x-ray energy spectral 

30 filter 42, and a high-energy two-dimensional detector 44. After the filter 42, x-ray energy is higher than that before the 
filter 42. The image detected by the low-energy detector 40 is denoted by D^rjJ) and the image detected by the high- 
energy detector 44 is denoted by D^fl j). Preferably, the low-energy x-rays have an average energy of from 10 keV to 
1 00 ke V and high-energy x-rays have an average energy of from 30 keV to 500 keV, with the high-energy x-rays having 
a higher energy than the low-energy x-rays. 

35 The purpose of measuring a pair of dual-energy primary images D^OJ) and D^iJ) of the rear detector assembly 

is to provide data for uniquely determining the corresponding primary low resolution image Dfp ( (x(i),y(j)) at the front 
detector. 

A flow diagram describing the method for determining a scatter-free image using the hardware of the first embod- 
iment is shown in Fig. 6. The x-ray source 14 emits x-rays with a uniform angular distribution and with an energy 
40 spectrum of <D 0 (E). The x-rays passing through the subject 1 2 carry information on the thickness and material compo- 
sition of the subject 1 2, expressed as projection mass density in units of grams/centimeters 2 (g/cm 2 ). The image induced 
by the x-rays incident on the front detector 1 6 is denoted as D fh (x,y) and is 

45 Dfh(x,y) = J [<£> 0 (E) x exp<-(p, b (E) x b(x,y) + (2) 

(i. s (E) x s(x,y))] x S f (E)dS + 
J O s (E) x S f (E)dE 

so 

where |i b (E) is the mass absorption coefficient of bone tissue and u^(E) is the mass absorption coefficient of soft 
tissue, with both n b (E) and u^(E) expressed in units of centimeter 2 /gram (cm 2 /g). Both of these values are known, 
having been determined experimentally and tabulated many years ago. S^E) is the x-ray spectral sensitivity (the elec- 
trical signal amplitude from the detector as a function of the number of x-rays with energy E incident upon the detector) 
55 of the front detector 16. The term J 0> S (E) X S f (E)dE represents the signal caused by scatter. 

The rear detector assembly 26 has two detectors 40, 44, so there are two low-resolution images D^fiJ) and D^, 
(i,j), which are 
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D rL i(i,j) = J [<t> 0 (E) x exp(-(jx b (E) x b ( i , j ) + (3a) 
|J. S (E) x s<i,j))] x S rL (E)dE 



and 



DrHi(i,j) = J [<D 0 (E) x exp(-(ji b (E) x b ( i , j ) + (3b) 
|I S (E) x s(i,j))] x S rH (E)dE 



Equations 3a and 3b constitute a simultaneous equation system, where the values lor the signal pair D^i.j), D^, 
(i,j) are known quantities, being evaluated from the electrical signals of the rear detectors 40, 44. b(i,j) and s(i,j) are 
the unknown quantities for which equation pair 3a, 3b must be solved, as described below. Generally speaking, in 
mathematics, such a nonlinear simultaneous equation system may have an infinite number of solutions, may have 
multiple-value solutions, or may not have any solution. However, for the specific physical case of x-ray imaging, where 
the energy range is limited, preferably between 10 keV and 500 keV, it can be mathematically proved that a unique 
solution always exists for the nonlinear equation system 3a and 3b. 

Accurate b(i,j) and s(i,j) are calculated by the data decomposition method of the present invention, as will be 
described below. Now that the values for b(i,j) and s(i,j) are known, the front low-resolution scatter-free image D^^x, 
y) can be obtained for those front detector cells (x(i).y(j)) that are on the selected projection lines. Dfp|(x(i),y(j)) is the 
signal induced by only primary x-rays at the detector cell (x(i),y(j)) on the front detector, and is 

D f?1 (x(i) ,y ( j ) ) = J [<J> 0 (E) x exp(-(n fa (E) x (4) 

b(i, j ) + ji s ( E ) x s (i, j ) ) ] x 
S f (E)dE 



where (x(i),y(j)) is the coordinate of the front detector cell (x,y) lying on the same selected projection line as the rear 
detector cell (i,j) and the energy dependent function 0 O (E) x S^E) is given in calibration, as will be described in the 
data decomposition section below. 

Next, the low-resolution front scatter image D fS ,(x(i),y(j)) iS determined by applying equation 3, 



D fSI (x(i), (j)) = D fl (x(i),y(j)) - D fpi (x(i),y(j)) 

Because of the scatter is largely caused by Compton scattering, which has a rather smooth angular distribution, 
the low-resolution scatter image DfsrfxfQ.yQ)) can be extended to the entire (x,y) plane through interpolation without 
losing accuracy, yielding the high- resolution scatter image D^x.y). The high-resolution scatter image D|s h (x,y) is 
subtracted from the experimentally measured image D^fx.y), yielding the high-resolution scatter-free signal Dfp h (x,y). 

It is now clear why it is crucial to acquire a pair of dual-energy x-ray images D^i.j), D rH (i,j) to determine the scatter. 
The primary signals at the front detector Dfp|(x(i),y(j)) can be uniquely determined only when the material composition 
image pair b(i,j) and s(ij) is known (see equation 7). If only one image at the rear detector is used, it would always be 
found that the signal ratio between the front detector and the rear detector is dependent on the energy spectrum or 
dependent on the image subject. Consequently, the primary image on the front detector could not be determined. 

It is also noted that, as a supplement to this general rule, there are two special cases where the relationship 
between the signals on the image plane (i,j) and the image plane (x(i).y(j)) can be degenerated or simplified to be 
proportional to each other and independent of the image subject. Thus, for correcting scatter, it is not necessary to use 
the dual-energy method for these two special cases. 

The first special case occurs when the range of the x-ray energy spectral distribution is sufficiently narrow so that 
the x-rays can be approximated as having only a single energy or a well-defined average energy E 0 . The signals on 
the front detector and the signals on the rear detector become a constant .independent of image subject and can be 
predetermined before using the system for imaging operations. 
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The second special case is shown in Figs. 7a and 7b, where only one detector is used for acquiring both the high 
resolution image D fh (x,y) and the low resolution image D H (i,j). A mechanical device moves the beam selector to and 
away from the detector's front surface. When the beam selector is present, as in Fig. 7a, blocking the scatter from 
reaching the detector, a low resolution image D rt (iJ) is acquired. When the beam selector is absent, as in Fig. 7b, 

s allowing all x-rays to reach the detector, a high resolution image D fh (x,y) is acquired, a portion of which is D fh (x(i),y(j). 
Here the subscript f denotes that the beam selector 18 is in its transmissive position not in front of the detector 16 and 
the subscript r denotes that the beam selector 18 is in its blocking position in front of the detector 16. In this special 
case, when the rest conditions are maintained substantially identical, if, and only if, the x-ray beam transmission through 
the holes is 100%, an energy-independent constant exists, which is independent of the image subject 12 and is pre- 

10 determined when the image subject 12 is absent. 

D fPI (x(i),y(j)) = C g (i,j) x D rt (i,j) (5) 

Fig. 8 shows the data processing procedures for the second special case. 
Second embodiment 
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In this embodiment, shown in Fig. 9, the front detector assembly 16 and the rear detector assembly 26 each have 
only one detector The x-ray source 14 emits two consecutive pulses, a low-energy pulse followed by a high-energy 
pulse. 

As shown in Fig. 10, the acquired first image pair includes the high-resolution images from the front detector 16 
and are 

Bf H h(x,y) = J [<X>oh(E) x exp(-(|i b (E) x b(x,y) + (6a) 
fl s (E) x s(x,y))] x S f (E)dE + 

J O fs (E,x,y) x S f (E)dE ■ , 

and 

De L h(x,y) = J [<E>ol(E) x exp(-(Ji b (E) x b(x,y) + ;6b) 
fi s (E) x s(x,y))] x S f (E)dE + 
/ O fs (E, x, y) x S f (E) dE 

and the second image pair includes the low-resolution images from the rear detector 26 and are 

DrHi(i.j) = J [Ooh(E) x exp(-(ji b (E) x b ( i , j ) + (7a) 
l± s ( E ) x s ( i 7 j ) ) ] x S r (E) dE 

and 



DrLi(i,j) '= I [<&ol(2) x exp(-(^i b (E) x b ( i , j ) + (7b) 
55 M- S (E) X s(i,j))] x S r (E)dE 

In the equation pair 7a, 7b, the acquired low-resolution image data are free of scatter radiation. By using the dual- 
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energy data decomposition methods described below, the simultaneous equation pair 7a. 7b is solved to find the 
solutions for the image pair of material composition b(i,j) and s(i,j). 

As described above with reference to the first embodiment, because the rear detector cell (i,j) and front detector 
cell (x(i),y(j)) lie on the same selected projection line, the low-resolution front detector primary image pair D fHP ,(x(i),y 
(j)). DfLPiMO.yG)) can be determined from the rear detector primary image pair D^i.j), D^fi.j). The front detector 
scatter image pair D fHSI (x(i),y(j)), D fLS ,(x(i),y(j)) are found by the equations 



D fHS ,(x(i).y(j)) = D h (k(I). y(j» - D (HPI (x(i),y(j)) (8a) 

w 

and 



DfLSiWO.yG)) - D fU (x(i),y(D) - D^MO.yffl) (8b) 

15 

As above, the low-resolution scatter images can be extended to all front detector cells not on the selected projection 
lines through interpolation without loss of accuracy to yield the high-resolution scatter image pair D^s^x.y), D fLSh (x, 
y). The high-resolution scatter-free images on the front detector assembly are denoted as DfHp^x.y) and DfLp^x.y) 
and are 

20 

D tHPh( x > y) = D fHh( x *y) - D fHSh( x «y) < 9a ) 

and 

25 

D fLPh< x -y) = D fu,< x .y) - D n.sh( x «y) <»> 

3Q The image pair D^phfx.y), D fLPh (x,y) is a pair of dual-energy x-ray images without scatter. This image pair in turn 

relates to the material composition of the subject by the equations 

D £H ph(x,y) = J [O 0H (E) x exp(-(|i b (E) x b(x,y) + (10a) 
35 |l s (E> x s(x,y))] x S«{E)dE 

and 



DfLPhfx.y) = J [Ool(E) x exp(-(u b (E) x b(x,y) + (10b) 
|I S (E) x s(x,y))] x S f (E)dE 



Thus, in addition to providing one scatter-free image, this embodiment provides a pair of scatter-free dual-energy 
images in the equation pair 10a, 10b. This equatbn pair is the fundamental dual-energy x-ray imaging equation system 
with the unprecedented feature that scatter radiation has been essentially removed from the two-dimensional detector 
In the equation pair 10a, 10b, the values of D fLPh (x,y) and D^p^x.y) are known from the above described calculations 
conducted on the image pair D fHh (x,y) D fLh (x,y) directly measured from the front detectors 16, and on the image pair 
D ru(U). DrHiOJ) directly measured from the rear detector 26. The unknown values are the two material composition 
images b(x,y) and s(x,y). 

The dual-energy x-ray data decomposition method can be further applied to the equation pair 1 0a, 1 0b. As a result, 
by using the quantitative relationships b = b(D H .D L ) and s = s(D H ,D L ) provided by the data decomposition method, a 
pair of high-resolution images b(x,y) and s(x,y) are readily obtained point by point for all front detector cells (x,y). The 
solution of the two-component material composition images b(x,y) and s(x,y) has a spatial resolution as high as the 
front detector 16 can provide. 

An alternate to the second embodiment substitutes an x-ray source generates x-rays continuously, alternating 
between high-energy x-rays and low-energy x-rays. 
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Another alternate to the second embodiment Inserts an x-ray energy filter between the x-ray source and the subject 
at the moment the x-ray source switches to generate the high-energy x-rays. The synchronization between the insertion 
of the filter and the switching high-voltage or double-pulse is preferably implemented by using a motor drive. The filter 
absorbs more of the low-energy x-rays, resulting in an increase in the energy difference between the low-energy x- 
rays and the high-energy x-rays. 

This invention not only provides a method and apparatus for removing scattering from two-dimensional detectors, 
but at the same time also provides a method and apparatus for scatter-free dual-energy x-ray imaging using two- 
dimensional detectors. 

The interrelationship between the method for removing scatter radiation in two-dimensional detectors and the 
method for dual-energy x-ray imaging using two-dimensional detectors can be summarized as follows: 

1. The method for removing scatter radiation from two-dimensional detectors utilizes and hinges on the method 
of dual-energy x-ray imaging free of scatter. Without dual-energy x-ray imaging, the scatter radiation cannot be 
accurately removed. 

2. The method for dual-energy x-ray imaging using two-dimensional detectors utilizes and hinges on the method 
of removing scatter radiation. Without substantially removing scatter from two-dimensional detectors, the accuracy 
of dual-energy x-ray imaging would be so degraded as to be meaningless. This invention solves both problems in 
a unified system. 

Third Embodiment 

The third embodiment, uses a pair of two-dimensional detector assemblies. The x-ray source 14 is a constant 
potential x-ray source that emits steady state x-rays, single pulse x-rays, or repetitive pulse x-rays with the same energy 
spectrum. The front detector assembly 16 has a low-energy two-dimensional detector 50, an x-ray energy spectral 
filter 52, and a high-energy two-dimensional detector 54. The rear detector assembly 26 also has a low-energy two- 
dimensional detector 56, an x-ray energy spectral filter 58, and a high-energy two-dimensional detector 60. The filters 
52, 58 operate in the conventional manner as described above with reference to the first embodiment. The front high- 
energy detector 54 is sensitive to higher x-ray energies than the front low-energy detector 50 and the rear high-energy 
detector 60 is sensitive to higher x-ray energies than the rear low-energy detector 56. 

This embodiment requires only single-pulse, constant-energy x-rays, as in the first embodiment, rather than the 
dual-energy x-ray source of the second embodiment. As shown in Fig. 12, following a single-pulse x-ray, two pairs of 
dual-energy x-ray images are acquired. The first pair includes the high-resolution images Dfl_ h (x,y) from the low-energy 
front detector 50 and D{H h (x,y) from the high-energy front detector 54. 

The second pair includes the low- resolution images D rU (i,j) from the low-energy rear detector 56 and D rH ,(i,j) from 
the high-energy rear detector 60. 

As described above with relation to the second embodiment, the image pair D^i.j) and D^i j) does not contain 
scatter. Thus, a pair of material composition images b(i,j), s(i,j) can be calculated for each cell in the (i,j) plane. This 
image pair is used to first determine the low-resolution primary image pair D fHP ,(x(i),y(j)), D fLPI (x(i),y(j)) and then to 
determine scatter image pair D fHS ,(x(i),y(j)), D ILSI (x(i),y(j)) in the same way as for the equation pair 14a, 14b. Scatter 
image pair D ms ,(x(i),y(j)), D fLS ,(x(i) J y(j)) is used to further determine the high -resolution scatter image pair DfH^x.y), 
D fLSh( x ' v )> as described above in the second embodiment. By subtracting the calculated high-resolution scatter image 
pair D fHSh (x,y), D^fx.y) from the front detector 16 high-resolution image pair D^fx.y), D fLh (x,y), the high spatial 
resolution image pair free of scatter D fHPh (x,y) and D fLPh (x,y) is numerically obtained. Using the inversion method 
described below, a pair of high-accuracy, high-resolution material composition images b(x,y) and s(x,y) can be further 
obtained. 

Data Decomposition Method 

The following is a step-by-step description of the data decomposition method summarized above and shown in 
Fig. 3. 

The first step is to construct the two simultaneous numerical equations D L = D L (b,s) and D H = D H (b,s) in three- 
dimensional space. A preferred method for doing this is to determine the function S(E) through the standard absorption 
method. 

An absorption curve is measured by using a collimated narrow primary x-ray beam. The electrical signal from a 
single detector cell D(t) as a function of the absorption plate thickness t is experimentally determined and is related to 
S(E) through the equation 
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D(t)= J S(E)X <D 0 (E) x exp(-n(E) X t)dE (11) 

Since the mass absorption coefficient ji(E) of the absorptbn plate material is known, the function S(E) can be 
5 determined. This method is especially convenient for the internal conversion type of two-dimensional x-ray detectors. 
In these detectors, the detection efficiency and detector energy response function can be expressed in a simple ana- 
lytical expression with few unknown parameters to be solved, 

io S(E) = {S 0 [l-exp( -Ho (E) x d) ] x aE} x (12) 

exp{-p. 1 (E) x d L - jll 2 (E) x d 2 ) 
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where S 0 (E)= [1-exp(-p 0 (E) x d)] x aE is the electrical signal amplitude induced by x-ray photons with energy E, |i 0 
(E) is the mass absorption coefficient of the detector's conversion layer, d is the thickness of the conversion layer of 
the detector cell, and where S^E) = exp'-u-^E) X d-, - ^(E) x d 2 ) is the x-ray transmission after leaving the image 
subject to the detector surface, m(E) and u^E) are the absorption coefficients of two given materials, and d-, and d 2 
are the thickness values of these materials. 

When the x-ray energy spectrum <I> 0 (E) is separately measured, these unknown parameters a, d, 6 V and d 2 are 
determined by using standard least square parameter-fitting techniques through equation 12. Then, the function S(E) 
is obtained to a high degree of accuracy for a single cell. 

Once the value for S(E) is determined to the desired accuracy, the dual-energy signals as a function of the materia! 
composition of the subject can be calculated through the equations 

D H = J S(E) x <D 0 (E)x exp(-(ji b (E) x b + (13a) 
(l s (E) x s)dE 



30 and 



35 



D L = J s(E) x <J> 0 (E)x exp(-(M*>(E) x b + (13b) 
fi s (E) x s)dE 



where u. b (E) and u s (E) are the well-documented mass absorption coefficients for bone tissue and soft tissue, respec- 
tively. The mass suriace densities b and s are assigned values that sufficiently cover the real range of the subject 12. 

Another preferred method for constructing the quantitative explicit functions D L = D L (b,s) and D H = D H (b,s) is to 
conduct direct measurements of signals D L and D H at a number of selected b and s values. The entire functions Dl = 
D L (b,s) and D H = D H (b,s) are obtained from the directly measured data points by using standard two-dimensional 
interpolation algorithms. The interpolation in this case is valid because the functions D L = D L (b,s) and D H = D H (b,s) 
are continuous, smooth, and monotonous. 

The second step is to determine the material composition images b and s as functions of the image pair Dh, D l . 
The procedures for obtaining a simultaneous equation system for b(D H ,D L ) and s(D H ,D L ) are shown graphically in Figs. 
13a to 13d. To do so, the simultaneous equation pair 24a, 24b must be inverted. A preferred method of inversion is as 
follows: (1 ) as in Figs. 1 3a and 1 3b, assign a pair of values in the desired range to b and s corresponding to one of 
the coordinate points in the (b,s) plane so that b = b n , and s = s mJ where n = 0,1, 2.. .N, and m = 0,1 ,2... M. Typical N 
and M values are in the range of between approximately 50 and approximately 5,000. The larger N and M, the higher 
the accuracy of the results. However, the largest values for N and M are limited by the available capacity of computer 
memory and computing speed. From the two numerical equations representing the three-dimensional surfaces D H (b, 
s) and D L (b,s), determine a pair of D H and D L values so that D H [n,m] = D H (b = b n , s = s m ) and D L [n,m] = D L (b = b n , s 
= s m)» where D H [n,mJ and DJn.m] are two specific real numbers, and (2) as in Figs. 13c and 13d, replot the four 
numbers D H [n,m], DJn.m], b n , and s m to provide a pair of data points on the three-dimensional surfaces b(D H ,D L ) and 
s(D H ,D L ). The data point on the three-dimensional surface b(D H ,D L ) is D H = D H [n,m], D L = DJn.m], b = b n , and the 
data point on the three-dimensional surface s(D H ,D L ) is D H = D H [n,m], D L = DJn,m], s = s m . After going through all the 
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b = b n values (bo.b-,, b 2 ,...,b N ) and all the s = s m values (S 0 ,S 1> S 2 ,...,S M ) J the essential part of the inversion task is 
complete. However, for the purpose of storing the inverted arrays b = b(D H ,D L ) and s = S(D H ,D L ), the step sizes of D H 
= D H [n,m] and D L = DJn.m] must be adjusted. In the inverted space, D H and D L are basis coordinates. From the N x 
M data points, only J data points are selected for D H and only K data points are selected for D L , where J and K are 
approximately in the same range as N and M. In the final form after the second step, two two<Jimensional arrays are 
obtained and stored: b = b(D H ,D L ) and s = s (D^DJ, where D H = D H [j], D L = D L [k]; j=0,1,2,...,J, D H [j] > D H [j+1] and 
k=0,1,2,...,K, DJk] >DJk+1]. Two additional one-dimensional arrays D H [j]and DJk] are also stored. Arrays D H [j] and 
DJk] are saved so that accuracy as high as real number calculations can provide is maintained. 

The third step is to find the desired results from the input data according to the established equations. The desired 
values for b and s at each cell location is determined by inserting the available data pair (D^DJ into the numerical 
equations of step 2. Conversely, the desired values for D H , D L , or only one of them if only one is needed, at each 
discrete cell location is determined by inserting the available data pair (b,s) into the numerical equations of step 1. 

The final step is to maintain the accuracy of the values for b and s in order to maintain a continuous domain function. 
The following procedures ensure elimination of the errors in connection with finite steps in data processing. 

In step 1, in the process of constructing the equation pair for D H [n,m] = D H (b = b n , s = s m ) and DJn,m] = D L (b = 
b n> s = s m)» for eacn P air o f values of b n and s m , the D H [n,m] and DJn.m] are measured or calculated to an accuracy 
of real numbers. D H [n,m] and D L [n,m] are stored in computer as real number arrays. 

In step 2, the inversion process, including replotting in D H space and D L space, introduces no errors due to the 
data processing. The step sizes can be changed without losing any accuracy as long as values for Dh = D H [j] are 
selected that are exactly equal to one of the D H [n,m] values that satisfies the condition D H [j-1] > D H [j] > D H [j+1], and 
values for D L = DJk] are selected that are exactly equal to one of the DJn.m] values that satisfies the condition DJk- 
1] > DJk] > DJk+1]. 

In step 3, for each measured dual-energy signal data pair (D HEX , D LEX ), first to find out the cbsest j and k values 
according to the criteria: D H [j] > D HEX > D H [j+1] and DJk] > D LEX > DJk+1]. From the index values j and k, the closest 
b and s are first determined as b 0 = b 0 (D H [j], DJk]) and s 0 = s 0 (D H [j], DJk]). The following equations give b and s 
values to an accuracy as high as real number calculations can provide: 

• b =- b 0 (D H [ j] ,D L [k] ) +. t3b(D H/ D L ) /3D L ] DH=D Htjj; dl-duu X (14a) 

(D LEX - D L [k]) + [3b(D H ,D L ) /dDnlDH-DHCj]; DL=DLik] x 

(Dhex - D H [j]) + higher order terms 

and 

s = s 0 (D H [j] ,D L [k] ) + [ds (D H ,D L ) /3D L ] dh=D h[ j ) ; dl=dl(IcJ X (14b) 
(D LEX - D L [k]) + [3s <D H ,D L ) /3d h ] DH=OH{jl; DL=DL [ k] * 
(D H ex - D H [jJ) + higher order terms 

where the values for the higher order terms are found in standard calculus textbooks. 

Also in step 3, if the image pair D H and D L from a given material composition data pair (b ex ,s ex ) is to be found, D H 
and D L is obtained to an accuracy of real numbers by using similar standard Taylor expressions. 

Thus, the procedures described above provide methods for directly solving the nonlinear dual-energy x-ray imaging 
fundamental equation systems in its original form with arbitrary bremsstrahlung spectra at an accuracy as high as using 
real number analytical calculations can provide. 

The following is a list of possible further variations in the embodiments: 

(1 ) Some well-established computation tools, such as sorting algorithms or database procedures, can be used to 
carry out the inversion process described above. The use of any such software package to conduct the above- 
described inversion process is contemplated by the present invention. 

(2) All of the steps described above, including the data decomposition method and the scatter elimination method, 
can be combined together to various degrees, from combining any two steps to combining all the steps into one 
procedure. For example, image pair (D^D^) can be used directly to find (DfHp.Dip) without explicitly determining 
(b,s). One way of doing so is to construct a pair of quantitative relationships D^p = (D^.D^) and D^p = (D^.D^) 
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in a data base and storing them. Any such combination methods are contemplated by the present invention. 

It is intended that the scope of the invention be limited not by this detailed description, but rather by the claims 
appended hereto. 



1. A two-dimensional.x-ray imaging system for taking images of a subject, said system comprising: 

10 

(a) an x-ray source adapted to emit x-rays for passage through said subject, said x-rays including primary x- 
rays having their direction of travel unaltered by interaction with said subject and said x-rays including scatter 
x-rays having their direction of travel altered by interaction with said subject; 

(b) x-ray detection means for detecting said x-rays in two dimensions; 

75 (c) said subject being located between said x-ray source and said x-ray detection means; 

(d) said detection means including high-resolution detection means for detecting said primary x-rays and said 
scatter x-rays simultaneously; 

(e) beam selection means for blocking passage of substantially all of said scatter x-rays and permitting passage 
of a portion of said primary x-rays; 

20 (f) said detection means including low-resolution detection means for detecting only said x-rays passed by 

said beam selection means. 

2. The x-ray imaging system of claim 1 wherein said beam selection means is substantially composed of an x-ray- 
absorbent material having a plurality of holes, the axes of said holes being parallel to the direction of travel of said 

25 primary x-rays. 

3. The x-ray imaging system of claim 1 wherein said detection means includes a detector assembly and wherein said 
beam selection means has two positions, a transmissive position wherein said high-resolution detection means 
includes said detector assembly, and a blocking position wherein said low-resolution detection means includes 

30 said detector assembly. 

4. The x-ray imaging system of claim 1 wherein said high-resolution detection means includes a front detector as- 
sembly and said low-resolution detection means includes a rear detector assembly, and wherein said system 
comprises, in physical sequence from front to rear, said x-ray source, said front detector assembly, said beam 

35 selection means, and said rear detector assembly said subject being located between said x-ray source and said 

front detector assembly. 

5. The x-ray imaging system of claim 4 wherein said x-ray source is adapted to emit two x-ray pulses of different 



6. The x-ray imaging system of claim 4 wherein said rear detector assembly includes, in physical sequence from 
front to rear, a rear low-energy detector, a rear x-ray energy spectral filter, and a rear high-energy detector, said 
rear low-energy detector including a plurality of x-ray-sensitive detector cells arranged in a substantially rectangular 
matrix, said rear high-energy detector including a plurality of x-ray-sensitive detector cells, the arrangement and 
quantity of said rear high-energy detector cells being substantially the same as the arrangement and quantity of 
said rear low-energy detector cells. 

7. The x-ray imaging system of claim 6 wherein said front detector assembly includes, in physical sequence from 
front to rear, a front low-energy detector, a front x-ray energy spectral filter, and a front high-energy detector, said 
front low-energy detector including a plurality of x-ray-sensitive detector cells arranged in a substantially rectan- 
gular matrix, said front high-energy detector including a plurality of x-ray-sensitive detector cells, the arrangement 
and quantity of said front high-energy detector cells being substantially the same as the arrangement and quantity 
of said front low-energy detector cells. 

8. A method for taking a two-dimensional x-ray image of a subject using a two-dimensional x-ray imaging system, 
said subject being composed substantially of two materials that interact differently with x-rays, said imaging system 
including an x-ray source, beam selection means, and a two-dimensional x-ray detector assembly having a plurality 
of detection locations, said subject being between said x-ray source and said detector assembly, said x-ray source 
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being adapted to emit x-rays for passage through said subject, said x-rays including primary x-rays having their 
direction of travel unaltered by interaction with said subject, said x-rays including scatter x-rays having their direc- 
tion of travel altered by interaction with said subject, said beam selection means having a blocking position wherein 
passage of substantially all of said scatter x-rays to said detector assembly is blocked and passage of said primary 
x-rays to said detector assembly is permitted, said beam selection means having a transmissive position wherein 
passage of substantially all of said primary x-rays and said scatter x-rays to said detector assembly is permitted, 
and selected detection locations being those of said detection locations receiving only said primary x-rays when 
said beam selection means is in said blocking position, said method comprising: 

(a) illuminating said subject with said x-rays; 

(b) moving said beam selection means to said transmissive position; 

(c) acquiring a raw high-resolution image from said detection locations and processing said image l^ to 
normalize it and to subtract dark signals, yielding a high -resolution image composed of said primary x- 
rays and said scatter x-rays; 

(d) producing, from said image D^, a low-resolution image D fl representing said selected detection locations; 

(e) moving said beam selection means to said blocking position; 

(f) acquiring a raw low-resolution image l ri from said selected detection locations and processing said image 
l ri to normalize it and to subtract dark signals, yielding a low-resolution image D ri composed of substantially 
only said primary x-rays; 

(g) calculating a low-resolution primary x-ray image Dfp, by multiplying said image D ri by a predetermined 
constant that is independent of said subject; 

(h) calculating a low-resolution scatter x-ray image Dfs, by subtracting said image Dfp| from said image D fl ; 

(i) calculating a high-resolution scatter x-ray image D,3 h by extending said low-resolution scatter x-ray image 
D fS , to all of said detection locations through interpolation; and 

(j) calculating a high-resolution primary x-ray image Dfp h by subtracting said high-resolution scatter x-ray image 
D fSh from said high-resolution image D m ; 

(k) whereby said image Dfp h is a high-resolution, two-dimensional primary x-ray image of said subject after 
said scatter x-rays have been substantially eliminated, said image Dfp h having a resolution substantially equal 
to the highest resolution available from said detector assembly. 

A method for taking a two-dimensional x-ray image of a subject using a two-dimensional x-ray imaging system, 
said subject being composed substantially of two materials that interact differently with x-rays, said imaging system 
including, in physical sequence from front to back, an x-ray source, a front two-dimensional x-ray detector assembly 
having a plurality of front detection locations, beam selection means, and a rear two-dimensional x-ray detector 
assembly having a plurality of selected rear detection locations, said subject being between said x-ray source and 
said front detector assembly, said x-ray source being adapted to emit x-rays for passage through said subject, 
said x-rays having an energy spectrum that can be approximated as having a single energy Eq, said x-rays including 
primary x-rays having their direction of travel unaltered by interaction with said subject and said x-rays including 
scatter x-rays having their direction of travel altered by interaction with said subject, said beam selection means 
blocking the passage of substantially all of said scatter x-rays to said selected rear detection locations and per- 
mitting the passage ol said primary x-rays to said selected rear detection locations, said selected rear detection 
locations receiving only primary x-rays, and selected front detection locations being those of said front detection 
locations intersected by x-ray projection lines extending from said x-ray source to said selected rear detection 
locations, said method comprising: 

(a) illuminating said subject with said x-rays; 

(b) acquiring a raw high -resolution image l m from said front detection locations and processing said image l m 
to normalize il and to subtract dark signals, yielding a high-resolution image D m composed of primary x-rays 
and scatter x-rays; 

(c) producing, from said image D m , a low-resolution image D fl representing said selected front detection loca- 
tions; 

(d) acquiring a raw low-resolution image l rt from said selected rear detection locations and processing said 
image l ri to normalize it and to subtract dark signals, yielding a low-resolution image D ri composed of substan- 
tially only said primary x-rays; 

(e) calculating a low-resolution primary x-ray image D^, at said front detector assembly by multiplying said 
image D r) by a predetermined constant that is independent of said subject; 

(f) calculating a low-resolution scatter x-ray image D^, at said front detector assembly by subtracting said 
image D fP , from said image D fl ; 
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(g) calculating a high-resolution scatter x-ray image D fSh by extending said low- resolution scatter x-ray image 
D fS , to all of said front detection locations through interpolation; and 

(h) calculating a high-resolution primary x-ray image Dfp h by subtracting said high-resolution scatter x-ray 
image Df$ h from said high-resolution image D m ; 

(i) whereby said image D fPh is a high-resolution, two-dimensional primary x-ray image of said subject after 
said scatter x-rays have been substantially eliminated, said image Dfp h having a resolution substantially equal 
to the highest resolution available from said front detector assembly. 

10. A method for taking a two-dimensional x-ray image of a subject using a two-dimensional x-ray imaging system, 
said subject being composed substantially of two materials that interact differently with x-rays, said imaging system 
including, in physical sequence from front to back, an x-ray sou rce, a front two-dimensional x-ray detector assembly 
having a plurality of front detection locations, beam selection means, and a rear two-dimensional x-ray detector 
assembly having a plurality of selected rear detection locations, said subject being between said x-ray source and 
said front detector assembly, said x-ray source being adapted to emit x-rays at two different average energy levels 
for passage through said subject, said x-rays including primary x-rays having their direction of travel unaltered by 
interaction with said subject and said x-rays including scatter x-rays having their direction of travel altered by 
interaction with said subject, said beam selection means blocking the passage of substantially all of said scatter 
x-rays to said selected rear detection locations and permitting the passage of said primary x-rays to said selected 
rear detection locations, said selected rear detection locations receiving only primary x-rays, and selected front 
detection locations being those of said front detection locations intersected by x-ray projection lines extending from 
said x-ray source to said selected rear detection locations, said method comprising: 

(a) illuminating said subject with x-rays of a first average energy level H; 

(b) acquiring a raw high-resolution image \^ from said front detection locations and processing said image 
l^ to normalize it and to subtract dark signals, yielding an image D mh which is composed of primary x-rays 
and scatter x-rays; 

(c) producing, from said image DfH h , a low- resolution image D mi representing said selected front detection 
locations; 

(d) acquiring a raw low-resolution image l^, from said selected rear detection locations and processing said 
image l^, to normalize it and to subtract dark signals, yielding an image composed of substantially only 
said primary x-rays; 

(e) illuminating said subject with x-rays of a second average energy level L; i : 

(f) acquiring a raw high -resolution image l^ from said front detection locations and processing said image 
IfLh to normalize it and to subtract dark signals, yielding an image D tfLh which is composed of primary. x-rays 
and scatter x-rays; 

(g) producing, from said image D^, a low- resolution image representing said selected front detection 
locations; 

(h) acquiring a raw low- resolution image l^, from said selected rear detection locations and processing said 
image l,y to normalize it and to subtract dark signals, yielding an image composed of substantially only 
said primary x-rays; 

(i) calculating a pair of low-resolution primary x-ray images D^p, and D^p, at said front detector assembly 
from a low-resolution dual-energy primary x-ray imaging pair composed of said image D rH , and said image D^,; 
(j) calculating a low-resolution scatter x-ray image D^g, at said front detector assembly by subtracting said 
image D^p, from said image D^, and calculating a low-resolution scatter x-ray image D^g, at said front detector 
assembly by subtracting said image D mp ^ from said image D rHI ; 

(k) calculating a high-resolution scatter x-ray image Dft_ Sh by extending said low- resolution scatter x-ray image 
D nsi to 3,1 of said front detection locations through interpolation and calculating a high-resolution scatter x-ray 
image by extending said low-resolution scatter x-ray image D^g, to all of said front detection locations 

through interpolation; and 

(I) calculating a high-resolution primary x-ray image D fHPh at said front detector assembly by subtracting said 
image D m sh from sa 'd image D mh and calculating a high-resolution primary x-ray image D^p^ at said front 
detector assembly by subtracting said image Dfus h from said image D^; 

(m) whereby said images D jHPh and O fLPh form a high-resolution, two-dimensional, dual-energy primary x-ray 
image pair D fHPh , D fLPh of said subject after said scatter x-rays have been substantially eliminated, said image 
pair having a resolution substantially equal to the highest resolution available from said front detector assembly 

11. A method for taking a two-dimensional x-ray image of a subject using a two-dimensional x-ray imaging system, 
said subject being composed substantially of two materials that interact differently with x-rays, said imaging system 
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including, in physical sequence from front to back, an x-ray source, a front two-dimensional x-ray detector assembly 
having a plurality of detection locations, beam selection means, and a rear detector assembly assembly, said 
subject being between said x-ray source and said front detector assembly, said x-ray source being adapted to emit 
x-rays for passage through said subject, said x-rays including primary x-rays having their directbn of travel unal- 
tered by interaction with said subject and said x-rays including scatter x-rays having their direction of travel altered 
by interaction with said subject, said rear detector assembly including a rear low-energy two-dimensional x-ray 
detector having a plurality of selected rear low-energy detection locations, a rear x-ray energy spectral filter, and 
a rear high-energy two-dimensional x-ray detector having a number of selected rear high-energy detection loca- 
tions substantially equal to said plurality of selected rear low-energy detection locations, said beam selection means 
blocking the passage of substantially all of said scatter x-rays to said selected rear low-energy detection locations 
and said selected rear high-energy detection locations and permitting the passage of said primary x-rays to said 
selected rear low-energy detection locations and said selected rear high-energy detection locations, said selected 
rear low-energy detection locations and said selected rear high-energy detection locations receiving only primary 
x-rays, and selected front detection locations being those of said front detection locations intersected by x-ray 
projection lines extending from said x-ray source to said selected rear low-energy detection locations and said 
selected rear high-energy detection locations, said method comprising: 

(a) illuminating said subject with said x-rays; 

(b) acquiring a raw high -resolution image l m from said front detection locations and processing said image l m 
to normalize it and to subtract dark signals, yielding a high -resolution image D^, which is composed of said 
primary x-rays and said scatter x-rays; 

(c) producing, from said image D^, a low- resolution image D fl representing said selected front detection loca- 
tions; 

(d) acquiring a raw low-resolution image l,^, from said selected rear high-energy detection locations and 
processing said image l^, to. normalize it and to subtract dark signals, yielding a low-resolution image D^, 
composed of substantially only said primary x-rays; 

(e) acquiring a raw low-resolution image l^ from said selected rear low-energy detection locations and 
processing said image l^, to normalize it and to subtract dark signals, yielding a low-resolution image D^, 
composed of substantially only said primary x-rays; 

(f ) calculating a low-resolution primary x-ray image Dfp, at said front detector assembly from a low-resolution 
dual-energy primary x-ray imaging pair composed of said image D rH) and said image D^,; 

(g) calculating a low-resolution.scatter x-ray image D^, at said front detector assembly by subtracting said 
low-resolution primary x-ray image Dfp, from said low-resolution image D fl ; 

(h) calculating a high-resolution scatter x-ray image D^ h by extending said low- resolution scatter x-ray image 
Dfs, to ail of said front detection locations through interpolation; and 

(i) calculating a high-resolution primary x-ray image D fPh by subtracting said high-resolution scatter x-ray image 
D fSh from said high-resolution image D^; 

(j) whereby said high -resolution primary x-ray image D fPh is a two-dimensional image of said subject after said 
scatter x-rays have been substantially eliminated, said image D fPh having a resolution substantially equal to 
the highest resolution available from said front detector assembly. 

12. A method for taking a two-dimensional x-ray image of a subject using a two-dimensional x-ray imaging system, 
said subject being composed substantially of two materials that interact differently with x-rays, said imaging system 
including, in physical sequence from front to back, an x-ray source, a front detector assembly, beam selection 
means, and a rear detector assembly, said subject being between said x-ray source and said front detector as- 
sembly, said x-ray source being adapted to emit x-rays for passage through said subject, said x-rays including 
primary x-rays having their direction of travel unaltered by interaction with said subject and said x-rays including 
scatter x-rays having their direction of travel altered by interaction with said subject, said front detector assembly 
including a front low-energy two-dimensional x-ray detector having a plurality of front bw-energy detection loca- 
tions, a front x-ray energy spectral filter, and a front high-energy two-dimensional x-ray detector having a number 
of front high -energy detection locations equal to said plurality of front low-energy detection locations, said rear 
detector assembly including a rear low-energy two-dimensional x-ray detector having a plurality of selected rear 
low-energy detection locations, a rear x-ray energy spectral filter, and a rear high -energy two-dimensional x-ray 
detector having a number of selected rear high-energy detection locations equal to said plurality of selected rear 
low-energy detection locations, said beam selection means blocking the passage of substantially ail of said scatter 
x-rays to said selected rear low-energy detection locations and said selected rear high -energy detection locations 
and permitting the passage of said primary x-rays to said selected rear low-energy detection locations and said 
selected rear high-energy detection locations, said selected rear low-energy detection locations and said selected 
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rear high-energy detection locations receiving only primary x-rays, selected front high-energy detection locations 
being those of said front high-energy detection locations intersected by x-ray projection lines extending from said 
x-ray source to said selected rear low-energy detection locations and said selected rear high-energy detection 
locations, selected front low-energy detection locations being those of said front low-energy detection locations 
intersected by said x-ray projection lines, said method comprising: 



(a) illuminating said subject with said x-rays; 

(b) acquiring a raw high-resolution image l,^ from said front low-energy detection locations and processing 
said image l^ to normalize it and to subtract dark signals, yielding a high-resolution image D^, which is 
composed of said primary x-rays and said scatter x-rays; 

(c) producing, from said image D^, a low- resolution image Dfy representing said selected front low-energy 
detection locations; 

(d) acquiring a raw high-resolution image l^ from said front high-energy detection locations and processing 
said image l^ to normalize it and to subtract dark signals, yielding a high-resolution image DfH h , which is 
composed of said primary x-rays and said scatter x-rays; 

(e) producing, from said image D^, a low- resolution image D^, representing said selected front high-energy 
detection locations; 

(f) acquiring a raw low-resolution image l^, from said selected rear high-energy detection locations and 
processing said image l^, to normalize it and to subtract dark signals, yielding a low-resolution image D^, 
composed of substantially only said primary x-rays; 

(g) acquiring a raw low-resolution image l^ from said selected rear low-energy detection locations and 
processing said image l^, to normalize it and to subtract dark signals, yielding a low-resolution image D^, 
composed of substantially only said primary x-rays; 

(h) calculating a pair of low-resolution primary x-ray images D^p, at said front high-energy detector and D^p, 
at said front low-energy detector from a low-resolution dual-energy primary x-ray imaging pair composed of 
said image D,^, and said image D^; 

(i) calculating a low-resolution scatter x-ray image D^gj at said front high-energy detector by subtracting said 
image D^p, from said image D mi and calculating a low-resolution scatter x-ray image D^s, at said front low- 
energy detector by subtracting said image D^p, from said image D^,; 

(j) calculating a high-resolution scatter x-ray image by extending said low-resolution scatter x-ray image 
DfHs, to all of said front high-energy detection locations through interpolation and calculating a high-resolution 
scatter x-ray image DfL Sh by extending said low-resolution scatter x-ray image Dfl_ s , to all of said front low- 
energy detection locations through interpolation; and 

(k) calculating a high-resolution primary x-ray image D^p^ by subtracting said image DfH Sh from said' image 
and calculating a high-resolution primary x-ray image 0 fLPh by subtracting said image DfL Sh from said 
image D^; : 
(I) whereby said images DfH Ph and D^ph form a high-resolution, two-dimensional, dual-energy primary x-ray 
image pair D mph) 0 fLPh of said subject after said scatter x-rays have been substantially eliminated, said image 
pair having a resolution substantially equal to the highest resolution available from said front detector assembly 



19 



EP 0 782 375 A1 




21 



# # 

EP 0 782 375 A1 



Normalization of detector 
cells in order to use a single 
detector ceil to represent 
all detector cells 



Calibration with 
experimental data 



h 



For any data pair 
<b,s), find (D H ,D L ) 



Constructing a nonlinear 
equation system D H = 0 H (b,s) 
and D L = D H (b,s) 



Numerical inversion of 
the nonlinear equation 
system for D H and D L 



v 



Maintain an accuracy 
as high as using real 
number analytical 
calculations can 
provide at each step 



For any data pair 
(D H ,D L ), find (b,s) 



H 



Reconstructing an 
equation system b = b(D H ,D L ) 
and s = s(D H ,D L ) 



FIG. 3 



22 



EP 0 782 375 A1 



Illuminate subject with x-rays 



1 


r 


Retrieve l rHI and l fLI from rear detector and 
process to normalize and subtract dark 
signals, yielding image pair D rH , and D lL) 




f 


Solve image pair D fHI and D fLI to | 
determine subject materials b and s 1 






Retrieve image l fh from front detector 
and process to normalize and subtract 
dark signals, yielding image D (h 




f 


Determine scatter c< 
at detector cells on 
by inserting val 
and solvir 


Dmponent D (SI of D |h 
th« projection lines 
ues for b and s 
>g for D fSI 



Illuminate subject with x-rays | 




r 


Retrieve image D m (x,y) from front detector | 




r 


Retrieve D fHI (i,j) and D rLI (i f j) from 1 
rear detector I 


i 


r 


Solve image pair D rH| (i t j) and D rLI (i,j) to | 
determine subject materials b and s I 




r 


Determine scatter-free image D fp| (x(i),y(j)) I 
by inserting values for b and s I 




t 



Determine scatter image D ls ,(x(i),y(j)) by 
subtracting scatter-free image D fpi {x(i),y(j)) 
from portion of D |h on projection lines 



Interpolate D fS| for 
not on projectio 
scatter ii 


front detector cells | 
n lines, yielding 1 

* a 9 e D tsh J 


* i 


r 


| Subtract D fSh from D 


fh to yield image D <Pn | 



FIG. 4 



Interpolate DfSI( 
detector cells not < 
yielding scatter 


x(i),y(j)) for front 
>n projection lines, 

ima 9 e D ish( x ^) 




r 


Subtract D |Sh (x,y) from D fh (x,y) to 

yield image D (Ph (x,y) j 



FIG. 6 



23 




EP 0 782 375 A1 




FIG. 5 



24 




25 



• m 



EP 0 782 375 A1 




26 




EP 0 782 375 A1 



Illuminate subject with x-rays f 




f 


Move beam selector means I 
away from subject and detector | 




f 


I Retrieve image D (h 


(x,y) from detector J 




r 


Move beam selector means 1 
between subject and detector | 


y 


f 


Retrieve image D rl (i,j) from detector | 






Determine the low-resolution primary 
image D fpi <x(i),y(j)) = D fl (i,j) X C g (i,j), 
where C g (i,j) is a predetermined constant 
independent of the image subject 




r 


Determine the low-resolution scatter image 
D fS| (x(i),y(j)) by subtracting scatter-free 
image D fpl (x<i),y(j)) from D, h (x(i),y<j)) 




f 


Interpolate D fSI {x(i),y(j)) for front 
detector cells not on projection lines, 
yielding scatter image D fSh (x,y) 




r 


Subtract D fSh (x,y) from D m (x,y) to 
yield image D fPh (x,y) 



FIG. 8 



IHur^ 



Retrieve^n^ 



Retrieve D, H! (i,j) from rear detector | 



Illuminate subject with x-rays of energy L I 

T 

|^etnev^ma^ 



£^^Retr^^ 



Solve image pair D rHI (i,j) and D rLI (i,j) to 
determine subject materials b and s 



Determine scatter-free image pair 
x(i),y(j)) and D fH(M (x(i),y(j)) by 
inserting values for b and s 



Determine scatter image D IHSI (x(i) t y(j)) by 
subtracting scatter-free images D (HPI (x(i),y(j)) 
from portion of D fHh on projection lines and 

determine scatter image D fLSI (x(i),y(j}) by 
subtracting scatter-free images D fLPI (x(i),y(j)) 

from portions of D ILh on projection lines 



Interpolate D IHSI (x(i),y(j)) for front detector 

cells not on projection lines, yielding 
scatter image D (HSh (x,y), and interpolate 
D u.si( x (')ty(j)) * or f ron * detector cells not on 
projection lines, yielding scatter 

lma 9 e °RSh( x »y) 



Subtract D IHSh (x,y) from D fHh (x,y) to yield 
image D fHPh (x,y) and Subtract D fLSh (x,y) from 

°fLh( x »y) to y |e,d ifna 9 e D ftPh< x .y) 



FIG. 10 



BNSDOCID: <EP 0782375A1 I > 



27 



• # 



EP 0 782 375 A1 




28 




BNSDOCID:<EP n7A**7RA1 I > 



29 



EP 0 782 375 A1 



| Illuminate subject with x-rays | 




r 


Retrieve images D fHh (x,y) and D fLh (x,y) § 
from front detector | 




r 


i 


Retrieve images D fMI (ij) and D fLI (i,j) 
from rear detector 




r 


i 


Solve image pair D rHJ (i,j) and D rU (i,j) to 
determine subject materials b and s 


i 


f 




Determine scatter-free image 
pair D, LP ,(x(i),y(j)) and D fHPI (x(i),y(j)) by 
inserting values for b and s 


I 


i 


f 




Determine scatter image D fHSI (x(i),y(j)) by 
subtracting scatter-free images D mpi (x(i),y(j)) 
from portion of D (Hh on projection lines and 

determine scatter image D fLSI (x(i),y(j)) by 
subtracting scatter-free images D api (x(i),y(j)) 

from portions of D fLh on projection lines 




r 




Interpolate D fHSI (x{i),y(j)) for front detector 

cells not on projection lines, yielding 
scatter image D, HSh (x,y), and interpolate 

D fLsiW).y(i)) ,or * ront detector cells not 
on projection lines, yielding 
scatter image D ILSh (x,y) 


i 




Subtract D fHSh (x,y) Irom D fMh (x f y) to yield 
image D IHPh (x,y) and subtract D 1LSh (x,y) from 

D iLh( x »y) to y ie|d ima 9 e D aPh( x »y) 



FIG. 12 



30 



EP 0 782 375 A1 




b=b. 



— b(D H sD H In i m] 9 D L sDJn,m]) 



s=s 



P L [n t m] 



D M [n,ml 



^ s(D H =D H [n,m] t D L =:D L [n,ml) 



FIG. 13c 



FIG. 13d 



BNSDOCID: <EP 078237SA1 I > 



31 



# 



EP 0 782 375 A1 



3 



European Patent 
Office 



EUROPEAN SEARCH REPORT 



Application Number 

EP 96 30 9554 



DOCUMENTS CONSIDERED TO BE RELEVANT 



Category 



Citation of document with indication, where appropriate, 
of relevant passages 



Relevant 
to claim 



CLASSIFICATION OF THE 
APPLICATION (lnt.CL6) 



A 
A 



EP O 218 923 A (GENERAL ELECTRIC COMPANY) 

* column 1, line 15 - column 6, line 7 * 

MEDICAL PHYSICS, 

vol. 17, no. 5, September 1990, NEW YORK 
US, 

pages 866-875, XPO0O17O780 

B.K. STEWART ET AL.: "SINGLE-EXPOSURE 

DUAL-ENERGY COMPUTED RADIOGRAPHY" 

* the whole document * 

EP 0 184 247 A (PHILIPS PATENTVERWALTUNG 
GMBH) 

* abstract * 

* page 3, line 26 - page 7, line 34 * 

EP 0 105 618 A (THE BOARD OF TRUSTEES OF 
THE LELAND STANFORD JUNIOR UNIVERSITY) 

* abstract * 

* page 1, line 19 - page 2, line 19 * 

US 4 829 552 A (R.J. ROSSI ET AL.) 

* abstract * 

* column 1, line 30 - column 3, line 8 * 

US 5 440 647 A (C.E. FLOYD ET AL.) 

* column 1, line 14 - line 58 * 

* column 4, line 16 - column 6, line 43 * 

US 4 389 729 A (J. STEIN) 

* column 1, line 18 - column 2, line 43 * 

* column 5, line 31 - column 6, line 36 * 



-/- 



The present search report has been drawn up for all claims 



1,2,8-12 
1,2,8-12 



HO5G1/60 
G21K1/02 
A61B6/02 



1,2,8 



1,2,8-12 



1,2,8 



1.8 



1,2,8 



TECHNICAL FIELDS 
SEARCHED (lnt.CI.fr) 



H05G 

G21K 
A61B 
H04N 
G01N 



THE HAGUE 



Due of cea«>l«fi«i of the ceartfc 

17 March 1997 



Horak, G 



CATEGORY OF CITED DOCUMENTS 



X : particularly relevant if taken alone 

Y : particularly relevant if combined with another 

document of the same c 
A : technological background 
O : non-written disclosure 
P : intermediate document 



T : theory or principle underlying the invention 
E : earlier patent document, but published on, or 

after toe filing date 
D : document cited in the application 
L>: document died for other reasons 

& : member of the same patent family, corresponding 



DMOrw^in. —era moootc* -• i _ 



32 



EP 0 782 375 A1 



European Patent 
Office 



EUROPEAN SEARCH REPORT 



Application Number 

EP 96 30 9554 



DOCUMENTS CONSIDERED TO BE RELEVANT 



Category 



Citation of document with indication, where appropriate, 
of relevant passages 



Relevant 



CLASSIFICATION OF THE 
APPLICATION (Int.CLn) 



JOURNAL OF APPLIED PHOTOGRAPHIC 
ENGINEERING, 

vol. 9, no. 6, December 1983, SPRINGFIELD 
US, 

pages 184-195, XP002027697 

M.J. YAFFE ET AL. : "SCATTERED RADIATION 

IN DIAGNOSTIC RADIOLOGY: MAGNITUDES, 

EFFECTS, AND METHODS OF REDUCTION" 

* the whole document * 



1,2,8 



TECHNICAL FIELDS 
SEARCHED (Int.Q.fr) 



The present search report has been drawn up for alt claims 



THE HAGUE 



Dale of coB«vlHkai oTIfei 



17 March 1997 



Horak, G 



CATEGORY OF CITED DOCUMENTS 

X : particularly relevant if tilt en alone 

V : particularly relevant if combined with another 

document of the same category 
A : technological background 
O : urn-written disclosure 
P : intermediate document 



T : theory or principle ooderlying the invention 
E : earlier patent document, but published on, or 

after the filing date 
D : document cited in the application 
L : document cited for other reasons 

A : member of the same patent family, corresponding 



33 



# 



This Page is Inserted by IFW Indexing and Scanning 
Operations and is not part of the Official Record 

BEST AVAILABLE IMAGES 

Defective images within this document are accurate representations of the original 
documents submitted by the applicant. 

Defects in the images include but are not limited to the items checked: 

□ BLACK BORDERS 

□ IMAGE CUT OFF AT TOP, BOTTOM OR SIDES 

□ FADED TEXT OR DRAWING 

□ BLURRED OR ILLEGIBLE TEXT OR DRAWING 

□ SKEWED/SLANTED IMAGES 

□ COLOR OR BLACK AND WHITE PHOTOGRAPHS 

□ GRAY SCALE DOCUMENTS 

□ LINES OR MARKS ON ORIGINAL DOCUMENT 

□ REFERENCE(S) OR EXHIBIT(S) SUBMITTED ARE POOR QUALITY 

□ OTHER: 

IMAGES ARE BEST AVAILABLE COPY. 
As rescanning these documents will not correct the image 
problems checked, please do not report these problems to 
the IFW Image Problem Mailbox. 



